Phantom thick brane in 5D bulk 
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I. INTRODUCTION 



A discovery of the accelerated expansion of the Universe [J, [3, [3, 0, 0] and a great deal of work attempting to model 
this phenomenon allows to suppose that existence of a special class of scalar fields - phantom fields - is possible in 
^ the Universe [b']. The phantom fields are some type of matter with the violated weak or null energy conditions. One 
{3jr)_ of the ways of introduction of such fields is consideration of theories of a scalar field with a negative sign before a 
kinetic term (ghost fields) 0, [1, In using the ghost fields, we have obtained regular wormhole, particle-Hke 
and brane-like solutions. Another way consists in a consideration of braneworld models of dark energy which also 
allow violation of the weak or null energy conditions [ll|, • 

The accelerated expansion of the Universe is being supported in phantom models by the phantom fields filling all 
[ the Universe. In this paper we suppose that the phantom field really fills the Universe, and we use this assumption for 
T— I • obtaining a thick-brane solution in a 5D spacetime with two interacting phantom scalar fields. Previously, one of us 
' has already obtained a similar regular solution for two usual interacting scalar fields [l3| . Here we want to compare 
i the solutions from [13] with a solution to be found below, and to clarify the difference between thick-brane solutions 
with phantom and usual scalar fields. 
OO Also we want to compare our results with the results from [ij] where thick branes with one usual scalar field are 

under consideration (see also earlier work on thick branes [HI). It was shown there that corresponding solutions have 
^ \ AdS asymptotics, the potential V{(p) must be an alternating function, and also some fine tuning condition should be 
satisfied. 



II. EQUATIONS AND SOLUTIONS 



We start from the Lagrangian 



(1) 



where capital Latin indices run over 0, 1, 2, 3, 5, e = +1 for the usual scalar fields, and e = —1 for the phantom ones; 
ip, X are two interacting scalar fields with the potential 

V{^, X) = ^(^' - mlf + ^{x'- ^^If + - K), (2) 
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where Vq is some constant. The corresponding energy-momentum tensor is 



(3) 



We consider a problem when all variables depend only on the extra coordinate r. Then metric for a 4D-brane 
embedded in an external five-dimensional spacetime can be written in the form 



Using this metric, the Einstein equations will be 
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where the first equation was obtained by summation of (*) and (J.) components. 
The equations for the scalar fields are 

^" + 4-^' = ^[2x^ + Ai{^^~ml)], 
a 

x" + 4-x' = x[2v^' + A,{x'-ml)]. 
a 

Let us examine the system of equations ©-([8]) with boundary conditions 

a(0) = ao, 

a'(0) = 0, 

^(0) = ipo, (^'(0) = 0, 

X(0) = Xo, X'(0) = 0. 



The conditions pTI)l - (fT^ and equation ^ give the following value of the constant Vq: 

TA ^1/2 2\2 I ^2 / 2 2\2 , 2 2 

Let us choose the following numeric parameters for numerical analysis: 

ao^ipo = 1, Xo = VcTe, Ai 0.1, A2 = 1.0. 



(7) 
(8) 



(9) 
(10) 

(11) 
(12) 



(13) 



(14) 



The model with e = -|-1 was considered in |13| . Here we examine a case when e = — 1, and compare obtained 
results with the results from [l^. At numerical analysis of the system ([S])-® one should solve a nonlinear problem 
for calculation of eigenvalues of the parameters mi and m2. This problem is being solved by the shooting method 
which is described in detail in [l^. The obtained results are presented in Figs.[l][3l The solutions for the case e = — 1 
are shown by solid lines, and for the case e = -1-1 - by dashed lines. The following values of the masses were found: 
mi « 1.93927 and m2 « 1.97696852 for the case e = -1, mi « 2.5220566669937 and m2 « 1.85622511328021 for 
e = +1. 

Let us estimate asymptotic behavior of the solutions. One can sec from ([6]) that asymptotically 



a(-i) 
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1^0 > 0, 



fc( + l) = \ - 



1 /A^ 



'm4-yo ) > 0, 



(15) 



(16) 



where a^o is an asymptotic value of the metric function, and subscripts (-1) and (-1-1) refer to cases with phantom and 
usual fields, respectively. Then, using these expressions and looking for a solution of equations (HI)-® in the form 



ipKimi-d(p, x~5x, 
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FIG. 1: The scalar fields profiles for the cases e = — 1 (solid lines) and e = +1 (dashed lines). 
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FIG. 2: The profiles of the functions = o'/a for the cases e = — 1 (solid line) and £ = +1 (dashed line). 
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we have the following asymptotic form for the scalar fields: 



<5<^(-i) 
SXi-i) 



ipoc exp 
Xooexp 
ifoo exp 
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<^X(+i) = XooCXp 



- 2k 
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(17) 
(18) 
(19) 
(20) 



where ipoo , Xoo are asymptotical values of the functions. One can see from the obtained expressions that asymptotically 
(5(y9(_|-i-) and 5x(±i) ~^ 0^ i-^- the solutions are finite. 
In Fig. [3] the profiles of the energy density 



e(r) = 

for the phantom and usual scalar fields are shown. 



(21) 



III. STABILITY ANALYSIS 



Let us now study the dynamical stability of the above solutions against linear perturbations. To do this we have to 
write field equations following from the Lagrangian ([T]) keeping time dependence of metric and scalar fields functions. 
Then (*) , {^) and ([ ) components of the Einstein equations will be: 
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where differentiation with respect to t and r is denoted by a dot and a prime correspondingly. 
Introducing a new metric function a via 



the component ([) takes the form 

3a' = e {-Iff' - xx') ■ 

The corresponding equations for the scalar fields are 

[^3 + 2aip] - ip" - Aa'ip' = -ip [2x^ + Ai {ip^ - ml)] , 
[X + 2a^] - x" - 4a'x' - -X [2p^ + A2 (x' - ml)] . 



(25) 

(26) 
(27) 

We perturb the solutions of the system (p5ll -(f27 l) by expanding the metric function and scalar fields functions to 
first order as follows 



a = aQ{r) + ai{r) cos Lut, 
— ipo{r) + pi{r) cos ujt, 
X = Xo{r) + Xi{r)cosujt, 



(28) 
(29) 
(30) 



where the index indicates the static background solutions of equations ([5])- ([8]), and the index 1 refers to perturbations. 
Then one has from (12511: 



-3 i^'oVi +X0X1) ■ 



Using the last expression, one can rewrite equations ((26|) . (f27l) as 
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For existence of stable solutions it is necessary to provide positiveness of an eigenvalue uP' . If there is a negative 



eigenvalue ui then the solution will be unstable since then Lpi , xi 



will grow exponentially. In order to make 



clear this question, we will search for numerical solutions of equations p2p . p3p with boundary conditions 

(pi(0) = l, ^'i(0) = 0, Xi(0)-0, x'i(0)=0. (34) 



For a case e 



-1 one can find that at 



a 1.6419 the following regular solutions for ^pi,Xi exist (see Fig. 0]). It 
indicates that the solutions with phantom scalar fields are stable for the case under consideration. But for a case 
e = 1 we could not find any regular solutions of equations ([32)) . ([33]). It seems that there is no stable solutions for 
usual scalar fields (at least for the model's parameters we used here). 



IV. DISCUSSION AND CONCLUSIONS 



In our opinion, the main distinction of the phantom thick brane from a brane created by usual scalar fields consists 
in asymptotical behavior of the metric function a(r): for the phantom thick brane a{r) k, Oooe'^'"^''' 00, and for 
the usual brane a{r) « aooe~'^<+^^'' 0. We think that it happens due to the following reasons: let a solution has 
the following form near r = 0: 

a(r) = ao + 04— H , (35) 

ip{r) = ^po + ^p2^—^ , (36) 

X{r) = Xo + X2y + ••• • (37) 
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Their substitution in ^ gives the foUowing result 

a4--|(^^+X^). (38) 

It imphes that at e < 0, a4 > and one has (at least near r = 0) a steadily increasing function which was confirmed 
by the numerical calculations presented in Fig. [51 At e > 0, 04 < and one has a decreasing function. 

Another interesting feature of the obtained solution is that the energy density of the phantom brane at the origin 
of coordinates r = has an absolute maximum. At the same time, the usual brane has a local minimum at r = 0. 
Physically, it implies that in the phantom case the brane looks like a hill on the background of phantom matter with 
constant energy density. In case of usual brane a maximum of energy is concentrated near the brane on both sides. 
In our opinion, it may results in different scenarios of trapping of 4D matter on the phantom and usual branes. 

At comparison of our solutions with the results from |14| one might note the following differences: (a) the metric 
function a{r) increases exponentially; (b) the potential V{r) = V {(p{r) , x{i^)) does not change a sign (see Fig. [5]); (c) in 
consequence of the condition (b), there is no a fine tuning condition. 

In Ref. [13] such fine tuning condition is necessary for obtaining of a thick brane supported by one scalar field. This 
fine tuning condition consists in the following: using the metric 

dsl = e^'^^'-^Tj^.dx^'dx-' - e^P^'Uz^ (39) 

with the harmonic coordinate z, such that ^/gg^'^ = — 1, one should has the following fine tuning condition 

V'(oo) = 0, (40) 

where 



V{z)^ J ^V{zi)dz,^ J e«^(^iV(^i)dzi. (41) 
In our case the corresponding condition could be obtained as follows: from equations (O, ([6]) one has 

r 

a'{r) = -e [ a{r)v(ip{r),x{r)) dr, (42) 
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FIG. 5: The profile of the potential energy eV{ip, x) for the cases e — ~1 (solid line) and e = +1 (dashed line). 



and because of a'(oo) ^ 0, we have 

oo 

J a{r)v(^ip{r),x{r))dr^O. (43) 



This indicates that the nontrivial potential V{ip, x) does not change a sign as it can be seen from Fig. [S] 

So we have obtained the solutions describing a thick phantom brane in 5D bulk. There arc two main distinctive 
features of these solutions: a) they have different asymptotic behavior; b) the solutions with phantom scalar fields are 
stable, and the solutions with usual scalar fields are not (at least for the parameters Ai, A2, (p{0), x(0) we used here). 
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